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Abstract 

We study the Newton-like problem of minimal resistance for a two- 
dimensional body moving with constant velocity in a homogeneous rar- 
efied medium of moving particles. The distribution of the particles over 
velocities is centrally symmetric. The problem is solved analytically; the 
minimizers are shown to be of four different types. Numerical results are 
obtained for the physically significant case of gaussian circular distribu- 
tion of velocities, which corresponds to a homogeneous ideal gas of positive 
temperature. 

Keywords. Two-dimensional Newton-type problems of minimal resistance, 
temperature motion, gaussian distribution of velocities, calculus of variations, 
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1 Introduction 

In 1686, in his Principia Mathematica, Newton considered one of the oldest 
problems of optimal control. The problem consists of finding the shape of a 
body, moving with constant velocity in a homogeneous medium consisting of 
infinitesimal particles, such that the total resistance of the medium to the body 
would be minimal. Newton assumed that the collisions of the particles with 
the body are absolutely elastic, the medium is very rare, so that the particles 
do not mutually interact, and that the particles are immovable, i.e., there is 
no temperature motion of particles. Problems of this kind may appear in con- 
struction of high-speed and high- altitude flying vehicles, such as missiles and 
artificial satellites. 
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Newton solved this minimization problem in the class of convex axially sym- 
metric bodies with the axis parallel to the velocity of the body, of fixed length 
along the axis and fixed maximal cross section orthogonal to the axis. Convexity 
implies that each particle hits the body at most once, and this fact allows one 
to write down the explicit formula for resistance. An account of the problem 
originally considered by Newton can be found, e.g., in the book p. 

Since the early days of the calculus of variations, many modifications of 
Newton's problem of minimal resistance have been studied in the literature 
In 1990th the interest to Newton's problem revived. Interesting results 
were obtained when dealing with the minimization problem in various classes 
of bodies obtained by withdrawing or relaxing the conditions initially imposed 
by Newton: axial symmetry [2 13 El; convexity (still maintaining the single 
shock assumption) single shock condition |10l 

On the other hand, the problem was studied under the more realistic assump- 
tions of presence of friction at the moment of collision (so that the collisions are 
not absolutely elastic) [7], and of mutual interaction of particles [2]. In the 
present paper, the case of temperature noise of particles is considered. First, we 
obtain general formulas in the d-dimensional case, d > 2, and then study the 
case d = 2 in more detail. 

We suppose that a convex and axisymmetric body is moving in M.^ with a 
constant velocity V, in a homogeneous medium of moving particles; the distri- 
bution of particles over velocities is the same at every point. In fact, it is more 
convenient to assume that the body is immovable, and there is a flux of particles 
falling upon it. This picture will be taken in the sequel. The length h of the 
body along the axis is fixed, and the maximal cross section by a hyperplane 
orthogonal to the axis is supposed to be a (d — l)-dimensional ball of radius 1. 

In section |21 the preliminary analysis of the problem in the d-dimensional 
case is made. In section O the minimization problem for c? = 2 is solved. In 
section 0] we consider the physically relevant special case of gaussian circu- 
lar distribution of velocities. Analytical formulas for resistance are given, and 
results of numerical simulations are presented. 

2 Calculation of Pressure and Resistance in the 
General Case 

Consider a flux of infinitesimal particles in M'', d > 2. The density of flux 
and the distribution of particles over velocities are the same at each point, the 
distribution being given by a density function p{v). The pressure of the flux at 
a regular point x € dB of the boundary of a convex body B cM.'^ equals n{nx), 
where Ux is the outer normal to dB at x, 




(1) 
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here x_ :— min{a;, 0} and (-I-) means scalar product; and resistance of the body 
to the flux equals 

TT{n,)dn''-\x), 



R{B) 



(2) 



where 'H'''~^ means [d — l)-dimensional Hausdorff measure. The formulas 1^ 
and Q are obvious modifications of the corresponding formulas from 3j. 
In what follows, we shall assume that 

(i) the function p is spherically symmetric with the center —Vcd, where 
V > and is the dth coordinate vector, i.e. p{v) = a{\v + Vedl); besides 

(ii) the function a is continuous and monotone decreasing, and 



r a{r) r dr < oo. 



The body B is supposed to be convex, open, and symmetric with respect 
to the dth coordinate axis. By translation along this axis and by subsequent 
homothety, it can be reduced to the form 



B = {{x\xd) : \x'\ < 1, < x^ < , 



(3) 



B C M'', where a;' = {xi, . . . , Xd-i), /+ and /_ are convex negative non- 
decreasing functions defined on [0, 1). We shall imagine that the dth coordinate 
axis is directed vertically upwards. The height of the body is 



-/+(0)-/-(0). 



(4) 



We shall suppose that h is constant. 

At a point x+ = {x\ — /+(|x'|)) of the upper part of the boundary dB, the 
outer normal equals 



1 



(5) 



and due to the property (ii) of p, from |^ one finds that pressure of the flux at 
this point equals Tr(nx^) — —p+ (/+(|a;'|)) • where 



0,...,0, 



1 



(6) 



Similarly, the outer normal to dB at a point x^ — {x\ of the lower 

part of dB equals 7r(n^_) = p- (/i(|a;'|)) • n^_, where 



P-{u) 



0, 



From ijHJl, O, and (QJ one obtains 

{viU ± Vd)J 



0, - 



1 



P±[u) = ± / 



p{v) dv. 



(7) 



(8) 
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Using the assumptions (i) and (ii), one can show that the functions p+ and 
p_ satisfy the foUowing properties: 

(a) p± eCi[0, +00); 

(b) there exist P±{u)] 

(c) p±(0) = hm„^+oo p'±{u) = 0; 

(d) p'+{u)<p'_{u)<Q. 

(e) There exist u± > such that is strictly monotone decreasing on 
[0, u±], and strictly monotone increasing on [u±, +00). 

The proof of these properties is not very difficult, but rather technical and 
long, and will be presented elsewhere. 

Let us calculate the resistance R{B), using the formula The integral 
in the right hand side of ^ is the sum of two integrals corresponding to the 
upper and the lower parts of dB. Change the variable in each of these integrals, 
substituting dW^-^ix) for ^f'+{\x'\Y + Ida;' and for fL{\x'\Y + ldx\ and 
substituting tt^Hx) for T:{nx:^) and for T:{nx_), respectively. One obtains 

J\x'\<i \ F I / 

P-(/:(k'l))-f-/i(k'l)^, 1) dx\ 

x-\<i \ Fly 

and as a result of integrating, using the fact that the functions p±{f±{\x'\)) and 
±/j_(|a;'|) jpj are symmetric and antisymmetric, respectively, with respect to x\ 
one gets 

P+{f'+{\x'\))dx^+ f p_(/:(|x'|))dx'l ed. 

lx'|<l J\x'\<l J 




R{B) = aa-i {n+{f+) +n-{f-)) ■ ea, (9) 
where ad-i is the volume of a unit ball in K*^^^, and 

7^±(/)= f\±{f'{t))dt''-\ (10) 

JQ 

From the formulas J^, H1U|) . and 10} it is seen that the problem of finding 
the functions /+, /_, minimizing R{B), can be solved in two steps. First, find 
the values 

inf 7^±(/) 

feM{h±) 

in the class A4{h±) of negative convex functions / such that /(G) = —h±. 
Second, minimize the sum inf_A4(/i^) 7?.+ (/)+inf^(/j_) TZ^{f) provided that + 
/i_ = h. 

Let us fix the sign "+" or "— ", and introduce shorthand notation p± = 
p, h± = h, TZ± = TZ. 

The following auxiliary lemma is a consequence of the Pontryagin Maximum 
Principle 
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Lemma 2.1. Let A > 0, fh E Ai{h), = 0, and let for any point t of 

differentiability of fh, u ~ fh{t) be a solution of the problem 

t"^"^ p(?i) + A u min . (11) 

Then fh is a unique solution of the minimization problem 

inf 7^(/), TZ{f)^ f\{f'{t))dt''-\ (12) 

f<iM(h) Jq 

Proof, (cf. 221) • Indeed, for any f £ M{h), f ^ fh one has 

+ A/'(i) > t'^-2p(/,;(i)) + A/,;(t). 

Moreover, this inequahty is strict on a set of t of positive measure. Integrating 
both sides of this inequality over t e [0, 1], one gets 



1 



d- 1 



p(/'(t))dt'^'i + A(/(l)-/(0)) 



> ^ ^ p(/;;w)rft'-' + A(/.(i)-A(o)), 

and using that /(I) < = and /(O) = /^(O) = — /i, one obtains that 

7^(/) > 7^(^). □ 



3 Solution of the Problem in the Two-Dimensional 

Case 

From now on, we shall assume that d = 2. 

Using the property (e) of the function p, it is easy to prove that the problem 

p{Q) -p{u) 

> max 

u 

has a unique solution; denote it by . 
Also, denote 

B m^p^ . 

There may appear three different cases: 

(a) As A > i?, the unique solution of ((TT|l is u = 0. 

(b) As A = i?, there are two solutions: m = and u ~ . 

(c) As A < B, the solution u is unique, besides u > and p'{u) = —A. 
Consider these cases separately. 

(a) X > B. The unique solution of Hll|) is u = 0, hence, according to lemma 
mH / = is the solution of (O for /i = 0, and n{f) = p{0). 
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(b) X = B. There are two solutions: u — and u — u^, hence a function 
fh G A4{h), whose derivative takes the values and u", is the solution of (|12() . 
By virtue of convexity of fh, there exists to e [0, 1] such that /^(t) = as 
t e [0, to], and /^(i) = u" as i e [to, 1]. Thus, 

f / as t e [0, to] 

■'''^^^ ^ \ -;i + iiO-(t-to) asie[io, 1]. ^ ' 

Taking into account that = 0, one concludes that h < vP and to = 1 — h/vP. 

Further, one has 

nh)= rp{o)dt+ [\{u")dt^p{o) + ^{p{u'>)-pm. 

Jo Jto " 

Using the definition of B, one gets 

n{h)=p{o)^Bh. 

(c) X < B. There is a unique solution u of pif) . hence the function 

fh{t)^~h + ut (14) 
solves (|12f) . where h = u. Here the resistance equals 

7^(^) = . 

Define the function 



p(0) - Bu if < It < M°, 
if w > 



Thus, one comes to 



Lemma 3.1. The solution fh of J^l'^) is given hv H'^) . if h < vP , and by ^ 

if h > vP ; moreover, TZ{fh) =p{h). 

Reverting to the subscripts "+" and "— " and using lemma im one concludes 
that the problem of finding 

inf (7^+(^J + 7^-(^_)) =:Rih) 

amounts to the problem 

min p;j(z), where phiz) ^ p+{z) + p-{h - z). (15) 

G<z<h 

The functions p±{u) are continuously differentiable on [0, +oo), and their 
derivatives are monotone increasing, hence p'^iz), < z < h is also monotone 
increasing. 
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Using the property (d), one concludes that B+ > hence there exists a 
unique value > u'^ such that = — Consider four cases: 

1) h< ul; 

2) u1<h< u,; 

3) u» < ft- < + u'L; 

4) /i > + it° . 

In the cases 1) and 2) one has p'i^{z) < p'i^{h) = p'_^_{h) + B- < as < z < /i, 
hence z = ft. is the solution of p5|) . Therefore, the optimal value of ft_ is zero, 
andA_(i) =0. 

1) ft < One has ft+ = ft < , hence fh^ is given by ((T^ . with io = 
1 — ft/u*].. The optimal body is a trapezium, with tangent of slope of its lateral 
sides equal to u^j.. The minimal resistance equals R(ft) — p+{0) — i?+ ft + p_(0). 

2) w+ < ft < u*. One has fh+{t) — —ft + ftt, hence the optimal body is an 
isosceles triangle. Here R(ft) = p+{h) +p_(0). 

In the cases 3) and 4) one has p+(ft) > —B- > — -B+ — hence 
ft > w° . Further, one has = p'+{h) — > 0; on the other hand, 

p'f^{u\) = p\{vP^) - p'_{h - u°^_) < B- - B+ < 0. It follows that the minimum 
of ph is achieved at an interior point of , ft] , so u'^ < ft+ < ft, and fh^ (t) = 
(t~l)h+. 

i) < h < u^, + vP_. Denoting ft = max{0, ft — one has ft < hence 
p'hih) = p'+{h) - PUh -h)< p'+{h) + i?- < , 

therefore the minimum of ph is reached at an interior point of [ft, ft]. Thus, 
< ft- < ft - ft < and fh_{t) = -ft- if i e [0, 1 - ft_/M°]; fh.{t) = 
-ft- + - 1 + ft_/M° ) if t e [1 - h^/u°_, 1]. The optimal body here is 
the union of a triangle and a trapezium turned over. The tangent of slope 
of lateral sides of the trapezium equals —vP_. The minimal resistance equals 
R(ft) +P-(0) - B_ (ft - u,). 

4) ft > u» + One has p'^(ft - ) = p+(ft - w° ) + B- > 0, hence the 
minimum of ph is reached at a point of ft — vP_). Thus, ft_ > u'?., and 
fh-{t) = th_. The optimal body is a union of two isosceles triangles with 
common base and of heights ft+ and ft- defined from the relations ft+ + ft- ~ ft, 
p^(ft+) = p'_{h ), ft+ > ft- > The minimal resistance here equals 
R(ft) =p+(ft+)+P-(ft-). 

4 Gaussian Distribution of Velocities — Exact 
Solutions 

Suppose that the density p is gaussian circular, with the mean —Vej, and vari- 
ance 1, i.e. 

pv{v)^^e-i^^+''^-^\ (16) 

Here, the value V is allowed to vary, so we shall denote the pressure functions by 
p±{u, V) instead oip±{u). Fixing the sign "+" and passing to polar coordinates 
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V = {—r simp, —r COS If) in the formula (jSJ, one obtains 

p+{u, V)^ j j p+ir, ^, V) rdrd^ , (17) 

where x+ := max.{x, 0}, and p+(r, (p, V) is the gaussian density written in 
the introduced polar coordinates, 

p+{r,^,V) ^ ^e-5('-'-2^'^™^^+^'). (18) 

ZTT 

Next, fixing the sign " — " and introducing polar coordinates in a slightly different 
way, V = (— rsiniy9,rcosiy9), one obtains 



7'^ (cos (p + u sin p)_^ ^ 



p_ (u,v) = - 1 1 — — V— n — ^ ^) ^'^^'^'^ • (19) 



Here p- (r, ]/) is the same density H16|l written in these coordinates, 

P-{r,ip,V) ^ _Le-^('''+2'^^™^^+^'). (20) 

ZTT 

Combining the formulas (|17() . (|18|l . (|19|l . and H20() . one comes to the more general 
expression 

cos sin\2>0 

Passing to the iterated integral and integrating over r, one obtains 

f T.N , e-^'/2 /- (cosy' + usinyj)^ 
P±(m,V^H± / T"^ — 2 1{±V cos p)d<p, 

TT J 1 + 

cos (p+u sin (f>Q 

where 

Changing the variable t — p — arcsin(u / Vl + u^), one finally comes to 

P±(u,'t^)=± / cos"^ T l{±z{T,u,V))dT, (23) 



where z(t, -u, y) = F (cost — m sin r) / vr+??_ 

Using one comes to the asymptotic formula 

z vl + u^ 
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V 



Figure 1: Height h versus velocity V of the body 



On the other hand, using if^ . in the hmit V — > oo one obtains the asymptotic 
formula for pressure 



P+{u,V) 



1 



(l + o(l)), p-{u,V) = o{V^)- 



here the main term coincides, up to the factor V"^ , with the pressure in Newton 
(i.e. zero-temperature) case. 

Numerical simulations were done using Maple. The results are presented on 
figures □ El and El 

Graphs of the three functions, shown on figure 1, separate the V-h plane into 
four regions, which correspond to the four possible solutions. The lower function 
tends to 1 as ^ oo. As = 0, the lower, the middle, and the upper functions 

take the values a, a, and 2a, respectively, where a — \J + VS) /2 w 1.27. 

Let R(V, h) be the minimal value of resistance at given values V and h\ define 
reduced resistance by R(y, h) = R(K h)/V'^. The graphs of R(y, h) versus h are 
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U- 

12 
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R{V,h)8-\ 
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V = 0.1 

V = 0.2 
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V = oo 



Figure 2: Reduced resistance R{V, h) versus height h of the body 
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V = 2 

V = 3 
y = 5 

V = oo 



Figure 3: Reduced resistance R(V, h) versus height h of the body 
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l-h/2 iih<l 

1/(1 + iih>l. 



shown on figures |21 and |31 for different values of V . As F ^ oo, R(V, h) tends to 
R(cxD, h) = 

As y — > 0, R(V, h) goes to infinity, besides 
^27^- lim {VK{V,h)) = 



2 - h/a^ if h < 2a 
4/(4 + if/i>2a. 



5 Conclusions and Final Comments 

In this paper, we treat one of the earhest problems in optimal control: Newton's 
problem of minimal resistance. We have obtained a full analytical description 
of the case when a two-dimensional body moves through a rarefied medium of 
infinitesimal particles, whose velocities are distributed according to the gaussian 
law. From the physical viewpoint, such a medium is just a homogeneous gas 
of positive temperature, while the case of immovable particles is related to a 
zero-temperature gas. 

The analytical formulas of the two-dimensional problem considered in this 
work are quite complicated, even for Maple, but using the current computational 
power, numerical simulations are made within a few days of calculations. The 
study of the tree-dimensional case would be a natural next step. However, 
calculation of the resistance functions and the whole analysis are more involved. 
This question is under development, and will be addressed elsewhere. 
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